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In ocean general circulation models, such as Bryan-Cox 1 and other 
derived models, the dissipation term is usually taken to have two axes of 
symmetry, with eddy coefficients A L and A„ in the symmetry plane and 
symmetry axis direction respectively. A L is much greater than A H , and 
both are much greater than the molecular value. If we further impose the 

condition of complete isotropy, the viscous stress term must return to the 

usual Laplacian of the velocity multiplied by a scalar, in whatever 
coordinates the equations are written and in particular in the earth 
surface coordinates in which the codes are written. However, the viscous 
stress terms used by nearly all the models, which were taken from 
Kamenkovich, do not do this. The reason for this flaw is the 
Kamenkovich's linearization of the gauge matrix occurred too early, i.e., 
the substitution of the earth radius a for the radius r = a + z was made 
too early, so that some of the z-derivati ves that should appear failed to 
appear. 

In this paper, the correct form of the viscous terms are presented. 
Indeed, the practical consequences of the error is probably not too 

serious, since the omitted z-derivati ves in Bryan-Cox 1 are multiplied by 
A h , and are therefore smaller compared to the terms multiplied by A L . 

Nevertheless, we present this paper in the interest of consistency arid 
possible future use. In the following, we first show our results in 
rectangular coordinates; then a revised form for the turbulent viscous 
term in earth-surface coordinates will be derived. The detailed formulae 
are given in the appendixes. 


1. Molecular and Eddy Viscosity 
in Rectangular Coordinates 


The eddy viscous force for incompressible fluids is the divergence of 
the Reynolds' stress (Appendix I) 
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an assumption widely accepted today. In (1.1), subscripts i and j stand 
for the three Cartesian components. 
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where u, v, and w are the x, y and z components of velocity; A, A„ and A, 
are the turbulence viscosity coefficients and they can be functions of x L 
y, z. 


For 3-dimensional isotropic fluids, i.e. A = A H = A, = constant 
the Reynolds' stress becomes: 
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and the eddy viscous term is then 
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For incompressible fluids, 
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(1.4) becomes 
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which of course can also be obtained by directly inserting a constant 
turbulence coefficient into (1.2). 


Obviously, under the complete isotropic conditions, this eddy viscosity 
term has the same form as molecular viscosity term - Laplacian of the 
velocity, the only difference being that the eddy viscosity coefficient is 
much larger than the molecular value. 
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2. Eddy and Molecular Viscosity Terms 
in Earth Surface Coordinates 


In rectangular coordinates, we have seen that when the transve 
symmetric condition is replaced by the complete isotropic condition, 
formulae for the eddy viscosity term returns to its spherical symmetric 
form This is easily seen in rectangular coordinates because the gauge 

laSx is constant. Choosing X, longitude ♦, 'tV^’th^eaJth s^?a e 
along the earth radius from the earth surface we get the earth ' surface 

coordinates. These are curvilinear coordinates in which the Reynolds 
stress (Appendix II) and the corresponding viscosity terms have more 
complicated forms; this fact, however, does not change the previou 
conclusion, since any physical phenomenon should not depend on coordinates 

chosen. 

In earth surface coordinates, the components of the gauge matrix are 


h A = (z + a) cos <t> 
h 0 = (z + a) 
h z = 1 

h = (z + a) 2 cos <f> (2.1) 


where a is the radius of the earth. 

The force F a , divergence of the Reynolds' tensor R a0 , is now 
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whose two horizontal components are; 
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Inserting the gauge matrix (2.1), the Reynolds' stress, 
of incompressible fluids 


the condition 
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and its related formulae shown in Appendix III into eq. (2.3) we obtain 
set of viscosity forces in the earth surface coordinate: 
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where A 2 stands for the 2-dimensional (X and 0) Laplacian operator, and A, 
A h and A l are the three eddy viscosity coefficients. 

If we let A = A h = A l = constant, we will obtain the viscous forces 
under the isotropic conditions: 
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These results are the same as that obtained from the direct transformation 
of the Laplacian from ordinary spherical coordinates to earth-surface 
coordinates. These terms, in fact, are identical to the molecular 
viscosity terms except the eddy viscous coefficients. 

Since the thin shell approximation, z « a, is used, z is usually taken 
to be zero in the denominators of eq. (2.7); in addition, u z is much 
smaller than u^ and u A , so it and its derivatives are omitted in these 
formulae as well. Then the correct approximation for the turbulent 
viscosity should be : 
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with the constant viscosity coefficients. 

If we now impose the isotropic condition, the viscous forces (2 9) 
become 
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equations in Ref. 2, we obtain instead 
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which are different from eq. (2.10); in other words, the equations in 
Ket. 2, under the isotropic condition, do not return to a spherical 
symmetrical form. This is why the equations in Ref. 2 are not a proper 

approximation for the eddy turbulent viscosity terms, although the 
last terms 3 
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in (2.8) 



Ua 


and 
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in (2.9) 


are very small compared with other terms in the equation (due to the very 
large radius of earth and small A H ), they should nevertheless be retained 
in a complete set of eddy viscous forces. 

The cause of this incompleteness is that the gauge matrix was 
linearized according to thin shell approximation too early, so the terms 
involving z-deri vati ves fail to appear when they should appear. 
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where q stands for the coordinates X, 0, z and the summation index 7 rolls 
over X, 0, z components. We have the following expressions for its 

components : 
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III. Two formulae related to the incompressible fluids condition (2 4) 
are very useful and are derived below: ' 
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